% Define approximate molar mass of components
molmass_yeastextract = 310; % g/L
molmass_tryptone = 71;

% Total carbon molar mass
molmass_LB = 10/molmass_yeastextract + 5/molmass_tryptone;
close all;
%% 37C LB
tempmat = monod_maxOD_LB_37C;
smat = size(tempmat);

conc1 = tempmat([1 2 3 4],1);
odmax1 = nanmean(tempmat([1 2 3 4],2:4),2);
odmax_err1 = nanstd(tempmat([1 2 3 4],2:4),0,2);

% Linear scale
figure;
errorbar(conc1/100*molmass_LB*1000,odmax1,odmax_err1, 'ko-');

set(gcf, 'position', [0 0 400 300])
set(gca, 'fontsize', 20);
xlabel('Concentration (mM)')
ylabel('Max OD600')
box off;
xlim([0 125]);

% Log scale
figure;
errorbar(log(conc1/100*molmass_LB*1000),log(odmax1),odmax_err1./odmax1, 'ko');
set(gcf, 'position', [0 0 400 300])
set(gca, 'fontsize', 20);
xlabel('Log concentration (mM)')
ylabel('Log max OD600')
box off;
